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Abstract 

After a rediscussion of the vanishing theorems for helicity amplitudes in unbroken gauge 
theories, we study the case of spontaneously broken gauge theories at high energy. The 
vanishing theorems generalize to a definite pattern of m/.E-suppression of the amplitudes 
that vanish in the massless case, where E is the energy scale of the process and m is the 
mass of the gauge vectors. We use only elementary arguments, and as an application we 
show how these methods can be employed to understand some aspects of the effective W 
approximation in the polarized case. 



1 Introduction 



Scattering amplitudes play a central role in the study of gauge theories, both from the theoretical 
and the phenomenological side. 

From the formal point of view, the properties of scattering amplitudes have often provided clues 
to hidden symmetries and unexpected dynamical structures in gauge theories. For example it was 
realized long ago that the tree- level amplitudes are effectively supersymmetric pQ [2] , so that their 
structure is significantly constrained by the fact that they obey supersymmetric Ward identities 
[3j[4j. By making use of these relations, it is easy to show that the tree- level amplitude for n 
(massless) gauge vector scattering vanishes unless it involves at least two helicities of each sign. 
Schematically, one says that the helicity configurations (± + ...+) and (± — ) vanish while 

the maximally-helicity-violating (MHV) nonvanishing amplitudes are of the type ( h..+) 

and (+ H — •••—), with all momenta ingoing. The tree-level MHV amplitude for n-gauge boson 
scattering can moreover be expressed by a remarkably simple formula, which was discovered by 
Parke and Taylor [5] and later proven by Berends and Giele [5]. It was also found that part of this 
simplicity extends to the loop level, at least for N = 4 super Yang-Mills theory [7], and there has 
been a resurgence of interest in uncovering new properties of scattering amplitudes after Witten 
[8] reformulated gauge theory as topological string theory in twistor space. Recent developments 
include a new type of recursion relation [9] [TD] and further generalizations [ll]|12j. 

At the phenomenological level, scattering amplitudes are critical to the prediction of cross 
sections at high-energy colliders, for processes within and beyond the Standard Model. The 
efficient evaluation of amplitudes involving many quarks and gluons is particularly important 
at the Large Hadron Collider (LHC), in which multi-jet final states are copiously produced, see 
[13] -[16] for modern reviews. 

The central issue at the LHC is however to test the electroweak sector, in which the gauge 
symmetry is spontaneously broken. A relevant question to ask is therefore how much of the 
simplicity of the massless case is left in the massive case, at least when the energy scale E at 
which the system is probed is much larger than the mass scale m of the gauge vectors. Since 
there is no discontinuity in taking the limit m/E — > for a gauge-invariant amplitude among 
given transverse-helicity states^, we expect the tree-level helicity amplitudes that vanish in the 
massless case to be in general nonvanishing, but suppressed by some power of m/E. In particular, 
the MHV amplitude in the massive case is in general expected to be of the type (+...+). The 
purpose of this study is to understand analytically the pattern of these suppressions, focussing on 
the tree-level amplitudes. 

In trying to extend the usual arguments it is clear that the supersymmetric one, although 
very elegant, is not straighforwardly generalizable to spontaneously broken gauge theories^]. An 
alternative way to understand why some tree-level amplitudes vanish in the massless case is to 
take advantage of the spinor formalism, in which the polarization vectors are expressed in terms 
of spinor wavefunctions. Although it is possible to generalize the spinor formalism to the massive 
casqj, and although this has proved to be a fruitful approach for numerical implementation, it 
is also clear that it is not suitable for an analytic understanding of the cancellations due to the 

Although the concept of helicity is not frame-independent in the massive case, as discussed below. 
2 For a generalization of the supersymmetric argument to the massive case see e.g. |17j . 
3 Seee.g. Q10I]- 
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complication of the various expressions. We take a different approach: since the spinor formalism is 
equivalent to working in axial gauge, we simply start from the general properties of the amplitudes 
in axial gauge. Instead of making use of the spinor representation of the polarization vectors, we 
thus directly use their defining properties. We show that in this way it is possible to study the 
(m/E) "-behaviour of the helicity amplitudes in a simple and direct way. Although the concept 
of helicity is frame-dependent, to define the "+" and "— " states amounts to choose a basis in the 
space of polarizations. We show that for a particular choice of the reference vector that defines the 
axial gauge, namely for light-like axial gauges, the MHV amplitudes are of the type ( — h •••+) and 
(H — ) while the amplitudes among equal helicities stay exactly zero. In this sense, light-like 
axial gauges seem to be the optimal choice in order to simplify the computations in the massive 
case. 

One may argue that such a study is not phenomenologically necessary since experimental data 
are compared to a given model through very efficient numerical codes that do not rely on the tree- 
level approximation and do not need to know whether some helicity amplitudes are suppressed, 
since they can compute 'brute force' and sum over all the states. Nevertheless, in order to extract 
from data the physically relevant parameters and couplings it can be very useful to have an 
analytic although approximate control of the relevant processes. In this sense the Effective W 
Approximation (EWA) originally discussed by Kunszt and Soper and recently reconsidered in 
[2Tj . although obsolete for the true data analysis, can be a useful instrument to understand under 
which conditions the sensitivity to the interesting physics related to the electroweak symmetry 
breaking is maximized. As an application of our study, we show that our methods can shed light 
on a few issues raised in [21]. 

This work is organized as follows. We recall some basic properties of the massless case in 
Section [2] and we consider the massive case in the axial gauge in Section |3j In Section H] we discuss 
the case of covariant gauge, in Section \5\ we make some contact with the EWA, and we conclude 
in Section [HI 

2 Massless case 

We start by recalling the most elementary argument, to our knowledge, that explains why the 
(± + ...+) and (± — ) amplitudes vanish in the case of unbroken gauge theories, see e.g. [13] . 
Here and throughout the paper, we take all momenta to be ingoing. We employ the axial gauge 
k^A^ = (see Appendix |A]), and we denote by e^(p,k) the polarization vector corresponding to 
a spin-1 particle with momentum p. The defining properties of such polarization vectors are: 



Notice that this choice of the polarizations corresponds to the usual one of the spinor formalism 
in which k is called the 'reference vector' (see e.g. [13]- [IS]), and that this gauge is singular when 
k ■ p = 0. Using either the spinor formalism, or setting m = in the arguments that are given in 



k ■ e ± (p, k) = p ■ e ± (p, k) = e ± (p, k) ■ e {p, k) = 
e ± (p,ky = e^(p,k) 
e ± (p, k)* ■e ± {p, k) = -1. 



(2.1) 
(2.2) 
(2.3) 
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the subsequent Section for the massive case, one can easily show the additional properties: 



where p, q, k are generic light-like momenta having nonzero scalar product with each other. Be- 
cause of gauge invariance, eq. (12.61) implies that one is free to choose a different reference mo- 
mentum for each polarization vector, and moreover one can choose different reference momenta in 
different gauge-invariant sub-amplitudes. As shown below, these properties are enough to prove 
that the (± + ...+) and (± — ) amplitudes vanish at tree level. 

Before going on, an observation is in order. In the following we will analyze the various 
amplitudes by counting the number of 4-vectors at our disposal at the numerator, that are the 
polaritazion vectors e M and the momenta p^ that appear in the Feynman rules for triple vertices. 
One has to notice, however, that the vector boson propagator M^ u (p, k) in axial gauge (1A.2I) has 
additional terms besides the metric tensor g^ u , and thus contractions like e^J\f^ u (p, k)e u involve in 
general also terms of the type e • p. A crucial point is to notice that all these terms are not present 
when the reference vector is light-like, which will be the case of our interest, and that in general: 



e^fa, k) (pi, k) Af% (pa, k)...M a ^{p ni k) cfo, k) = (-l) n e( gi , k) ■ e(g 2 , k) + O(-) . (2.7) 



where E is the typical size of the components of the various momenta and of the reference vector 
k. In this Section, the reference vector is always taken to be light-like. Of course the results do not 
depend on the reference vector in the case of unbroken gauge theories because of gauge invariance, 
but the proof is less transparent if k 2 ^ due to the additional terms in (12. 7J . 

• Case (+...+): Since there are only triple and quartic vertices, the number of vertices at 
tree level is < n — 2, where n is the number of external legs and the equality holds if the 
vertices are all with three legs. Moreover, each vertex brings at most one momentum 4- 
vector (if triple, no momentum if quartic). As a consequence, each term of the amplitude 
must involve at least one Lorentz contraction of two polarization vectors. Since the helicities 
are all positive, all such contractions vanish in axial gauge thanks to (12. 4p . 

• Case ( — h ...+): As above, each term of the amplitude involves at least one Lorentz con- 
traction of two polarization vectors. All the terms e + {pi, k) ■ e + (pj, k) vanish again provided 
that the reference vector k is the same for all the positive helicities. We are then left with 
terms of the type e + (pi, k) ■ e~(pi, ki), that vanish by (12.51) if we choose k — p\. Of course 
the reference vector of e~ must be a different 4-vector k\ ^ k, but this can always be done 
thanks to (I2.6p . as already observed. 

The amplitudes (—...—) and (H — ) also vanish in full analogy, while there is no suppression 
in general in the amplitudes involving at least two polarizations of each type. Other useful results 
concern the amplitudes involving vector bosons plus two massless scalar or fermion legs, which 
also vanish when all the vectors carry the same helicity. In the literature, this is usually shown 
using the supersymmetric Ward identities (see e.g. [13j[14j). We present here alternative and 
more elementary arguments that only make use of the properties of the polarization vectors in the 
axial gauge, along the same line as the discussion above. 



e + (p, k) 
e + (p, k) 



e + (q,k) = 
e-(k,q) = 



(2.4) 
(2.5) 
(2.6) 
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Case [cfxp* + ...+): In the case of an amplitude with n vectors and two scalar external legs, 
the number of vertices at tree level is < n, and the equality holds when all the vertices are 
triple. We thus need to consider only the amplitudes in which all gauge vertices are triple, 
since the others give a vanishing contribution because of (|2.4p . 

Let us call p\ the momentum of one scalar external leg, and let us consider the sub-piece 
of the amplitude corresponding to its first scalar-vector branching, with the emitted vector 
eventually ending in j < n vectors with momenta {k\, ...kj}. Taking all the reference vectors 
to be equal to p\, this subamplitude is of the form: 

j 

(2 Pl + J2 k iTA l ,{k l ,...k j ,e(k 1 ,p 1 ),...e(k j ,p 1 )) (2.8) 

i=i 

where A^ is the subamplitude attached to the vector line. We now show inductively that 
this expression vanishesj. For j = 1, this is true because A^ki, e(ki,pi)) = e M (fci,pi) which 
is orthogonal both to p\ and to k x by definition (12.11) . For j > 1, let us consider the first 
splitting of our vector line into two vector lines (we know that there are only triple vertices). 
To conclude it is sufficient to show that in all the terms in which (12. 8 p splits there is a 
factor in which the subamplitude attached to a vector line is again contracted only with 
the reference vector p\ or with the momentum along the line. This is sufficient because 
eventually we will get to an external leg and the contraction with the polarization vector 
will give zero as in the case j = 1. 

Call in fact p = — Yll=i ^ Q = Yll=i hi an d k = J2l=r+i ^» the momenta ingoing in the 
triple gauge vertex. The expression in (I2.8P is proportional to: 

(2 Pl - pY U ^ Pl) V(p, q, ky° A p {q) A a (k) (2.9) 

where we denoted by A p (q) and A p (k) the subamplitudes attached to the new vector lines, 
and (see Appendix [A} : 

KAp,pi) = -9v» + — Jf^ ViiVv ( 0) 

V(p,q,k) Vpa = g au (k-py + g»P(p- q y + g^( q -ky. (2.11) 

By expanding ( 12. 9 p using (I2.10p and (12. lip it is easy to see that in all the various terms either 
A p (q) or A a (k) is contracted either with the corresponding momentum (q or k, respectively) 
or with the reference vector p\. This is true apart from the term involving g pa , but this 
gives no contribution because it necessarily involves also the contraction of two polarization 
vectors. 

As a consequence, these amplitudes vanish at tree level. 



4 Notice that this can be shown more formally using BRST invariance, in particular the fact that a longitudinally 
polarized gauge boson is the BRST transformation of a ghost field. 
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• Case (ipip* + ...+): Consider first of all the fermion wavefunctions that appear in the 
amplitude. With massless fermions, there are only two independent solutions of the Dirac 
equation]! that can be written down as: 

</>±(p) = i(l±7 5 Mp) , Pu{ v ) = H±{V) = ®. (2.12) 

By definition these wavefunctions are eigenvectors of rotations along the axis specified by p, 
under which they transform as: 

V±(p) Rfy±(p) = e^h±(Rtp) = e T *'V±(p) • (2-13) 

At the same time, as shown in eq. (1A.14I) . the transformation properties of e ± (k,p) under 
the same rotation are: 

e^p) itfe^p) = e^e^EfanSp) = e ±l6 e ± (R%p) . (2.14) 

On the other hand, since ^(k^p)] = 0, the four spinors /l + '~(k,p)tfj + ,-(p) satisfy the 
Dirac equation so that each of them must be a linear combination of ip+(p) and ip_(p). But 
the above transformation properties under rotations along p can match only in two cases^l, 
and thus the other two must vanish: 

4 + (k,p)^(p) = ^ + (p)4 + (k,p) = 0, (2.15) 
/t(k,p)ip + (p) = ip_(p)4~(k,p) = 

For the amplitude with n vectors and two massless-fermion legs with momenta p± and P2, 
the number of vertices at tree level is again < n, and again the equality holds when all the 
gauge vertices are triple. Moreover because of (12. 4p we know that, when all the vectors have 
the same helicity, we can restrict ourselves to the terms in which all the gamma matrices 
are Lorentz-contracted with a polarization vector and never with a momentum 4-vector, 
otherwise a term e • e must be present. To conclude it is enough to notice that the helicity 
of the fermion line is necessarily conserved, that is, only the cases "0+(Pi) jl... /hp+ijp?) and 
"0-(Pi) 4— M , -(jp2) can be nonvanishing. Using (12.151) we then see that all the terms in the 
amplitude vanish if we choose as reference vector pi in the case of right-handed fermions, or 
P2 in the case of left-handed fermions. 

In conclusion, in this Section we reviewed known results about the vanishing of helicity ampli- 
tudes at tree level in an unbroken gauge theory. In doing so, we made use only of basic properties 
of the polarization vectors in axial gauge together with other elementary considerations. As is 
well known, these results greatly simplify the computation of tree-level scattering amplitudes in 
massless gauge theories, making the axial gauged a very convenient one for practical computations. 
Notice also that in the case (+...+), in order to see the cancellation, it is simply enough to employ 
an axial gauge while in the other cases it is crucial to choose properly the reference vector (s). As 
discussed in Section [31 this is the basic reason why in the massive case the (+...+) amplitude is 
suppressed by more powers of m/E than the others. 

5 Recall that there is no difference between u(p) and v(p) up to normalization conventions. This has to be the 
case since the projectors ^2u(p)u{p) and ^2v{p)v{p) are both equal to jp. 

6 One can further show that / ± (fc,p) esc ip±(p)ip±(k) + 'i/j T (k)ip T (p), from which (|2.15[) follows, see e.g. [T3] . 
7 Or the spinor formalism, that is equivalent to an axial gauge. 
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3 Massive case 



In order to understand the high-energy behaviour of the above helicity amplitudes in a broken 
gauge theory, the first step is to generalize f l2.4p - fl2.6p to the massive case. We consider the case 
in which the typical center-of-mass (CM) energy E of the process is much larger than the mass 
m of the vector boson, i.e. e e m/£ C 1. For simplicity, we consider the case in which the mass 
m is the same for all the gauge vectors. Generally speaking, a tree level amplitude involving n 
external legs has mass dimension 4 — n, so that we can write, working in the CM frame: 

A^-E^e*. (3.1) 

Since no gauge- invariant tree- level amplitude should become singular in the limit m — > 0, we 
expect t > 0. Notice that this is true in a "physical" gauge like the axial one, and thus it must 
be true in general for all the gauge-invariant subamplitudes. On the contrary non-gauge-invariant 
subamplitudes can grow faster with energy, and this behaviour cancels out only in the gauge- 
invariant results. As discussed in Section HJ this is what happens in covariant gauges in which 
the longitudinal polarization vectors scale as E/m introducing a "bad high-energy behaviour" 
that cancels away only in the gauge-invariant expressions. For this reason the axial gauge is the 
most suitable not only for an analytic understanding of the high-energy behaviour of the helicity 
amplitudes in the massive case, but also from the computational point of view since it does not 
involve large cancellations among subamplitudes. 

We also expect the polarized amplitudes that vanish in the massless case to be suppressed in 
(13. ip by some exponent t > 0, that we want to compute, or to stay equal to zero. In order to do 
that, as noted in [21], an useful selection rule can be obtained by noticing that the Lagrangian of 
the Goldstones s and gauge fields is invariant under the reparametrization: 

— )• , s — > — s , m — > —m . (3.2) 

As a consequence, the sign of m is not a physical observable and any amplitude among given 
external states must be either even or odd under m — > —m, since iterating (13. 2 p twice gives no 
change, and thus t must be an integer in (13. ip . Moreover we see from (IA.4|) that the longitudinal 
polarization vector changes its sign under ( I3.2p . while the transverse ones are unchanged, and 
thus we can transform an even amplitude into an odd one and vice-versa by changing a transverse 
external state into a longitudinal one. Since the amplitudes involving only transverse states are 
even under (I3.2p . this implies that the polarized amplitudes must be even or odd under m — > —m 
depending on whether they involve an even or an odd number of longitudinal states. 

Another useful argument is the one of recovering the massless case in the limit m — > 0. In fact 
in this limit the theory approaches an unbroken gauge theory interacting with massless scalars, 
for which the results of the Section [2] hold. In particular, we deduce that all the amplitudes among 
transverse states that vanish in the massless case must have t > in (13. ip . Moreover, when we take 
this limit in amplitudes involving longitudinal states, we can trade the longitudinally-polarized 
vectors for massless scalars charged under the gauge grourfl 

Notice however that, differently from the massless case, our polarizations e + (p, k) and e~(p, k) 
do not describe the helicity eigenstates in a Lorentz-invariant way since the helicity itself is not a 

8 This is also evident if we look at the explicit form of the longitudinal polarization vectors (|A.4|) . that makes 
the equivalence theorem an identity. 
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Lorentz-invariant concept. The physical states described by the polarization vectors e ± (p,k) are 
elements of a basis in the space of polarizations, and in the massive case this basis depends crucially 
on the vector k^. For this reason, in order for a polarized amplitude to have a physical meaning, 
it is now necessary that the reference vector is the same for all the polarizations. Analogously, as 
discussed in Section HJ the physical state described by e + (p, k) is in general not the same as that 
described by the usual 'plus' polarization in a covariant gauge e + (p), defined below. As we will 
see, in the massive case the basis of axial-gauge polarizations with light-like reference vector are 
a particularly convenient choice. 

To generalize (I2.4p - (l2.6p . we start from the projector on the transverse subspace. Using the full 
projector f]A.2|) and the explicit form of the 4- vector part of the longitudinal polarization (IA.4[) . 
we can write the transverse projector as: 

52<$(p,k)e$(p } ky = e+(p,k)e;(p,k) + e;{p,k)e+(p,k) (3.3) 

= -9^ ~ Z p2k 2 {P 2k n k » + k2 V^Vv - pk(p^K + ktf> v )) , 

with p 2 = m? and for any value of k 2 . To generalize (12.41) . we contract (13. 3p with e+(q,k) and 
square the two sides of the resulting equation. We obtain: 

2[ e +(p, k) ■ e +(g, k)\ [e-(p, k) • e +( g> k)] = ~k 2 ^^ 2 • 

By continuity with the case q — >■ p, there can not be suppression in e~(p, k) ■ e + (q, k). We conclude 
that: 

e+(p,k)-e + (q 7 k)=0(^) (3.4) 

where E ^> m is the typical size of the entries of p^ and k^ in the CM frame. Notice that the 
size of the entries of k 11 does not have any physical meaning, since the projector 0A.2p is invariant 
under a rescaling of the reference vector by an arbitrary factor. For notational convenience, we 
choose the leading nonzero entries of k^ to be of order E. 

To generalize (12. 5 p we do the same as above using e~(k, q), and we find: 

e + (p,k).e-(k } q) = O(^) (3.5) 
e + (p,k) ■ e~(k,p) = 0. (3.6) 

We also need to generalize f !2.6p . This can be done by writing down the explicit expression for the 
polarization vectors and expanding it for small p 2 /E 2 . The result is: 

4(P, k) - ej(p, q) = 0(1)| + 0(f^)f + • (3-7) 

A final observation before going on. As in the case m = 0, a crucial point in the various 
arguments for the cancellation or suppression of the amplitudes is the analysis of the possible con- 
tractions of the polarization vectors. To this end, note that massive gauge theories in axial gauge 
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provide an additional possible source of four momenta which can contract with the polarizations: 
the scalar- vector propagators Af^. Using eq. flA.2j) . it is immediate to show that: 



e«(p,k)Af^ q ,k) = 0^). (3.8) 

We are now ready to extend the discussion of Section [2] to the massive case. 

• Case (+ + ...+): Again the number of vertices at tree level is < n — 2, where n is the number 
of external legs. It is also still true that each vertex brings at most one four momentum, 
and thus each term of the amplitude must be proportional to the Lorentz contraction of 
two polarization vectors, keeping (12 .7p in mind. Since the helicities are all positive, all such 
contractions are suppressed at the level of 0{k 2 /E 2 ) by (13. 4p and (12. 7p . But in this case 
there is no condition on k 2 , in particular we can choose k 2 = 0. This argument shows that 
this kind of amplitudes vanishes at tree level also in the massive case, provided that one 
chooses a light-like reference vector. 

• Case ( — h ...+): For continuity with the case m = 0, we already know that this amplitude 
is suppressed by e t with t > 0, and from the selection rule discussed above we know that t 
is an even integer and thus t > 2. Through an explicit computation one can then see that 
the suppression is indeed no more than this, i.e. t = 2. 

An intuitive way to understand why t = 2 can be to recall that, in the proof of the vanishing 
theorem for the case ( — h ...+) in the previous Section, we had to make a precise choice of 
reference vectors. In particular the reference vectors for the positive-helicity polarizations 
had to be all equal to the momentum of the negative-helicity leg. Since this momentum is 
now time-like, we expect the suppression to be no more than t = 2 because of (12. 7p . 

• Case (0 + ...+): From (1A.4f) we see that, with k 2 = 0, we have e°(p,k) oc k and thus 
all the products e°(pi, k) ■ e + (pi, k) vanish. As recalled in the Appendix lAl the longitudinal 
polarization has also a fifth 'scalar' component, but when this is involved there must be either 
a contraction e + (pi, k) ■ e + (pj, k) or £^(pi, k)Af^ s (q, k) (q being some linear superposition of 
the Pj), both of which vanish for k 2 = 0. This shows that these amplitudes vanish at tree 
level also in the massive case, provided that one chooses a light-like reference vector. 

• Case (0 — h..+): As above, choosing k 2 = we have e°(p,k) cx k and all the e + (pi,k) ■ 
e + (pj,k) vanish. There is however no suppression in the contractions e~(p2,k) ■ e + (pi,k). 
In this case, taking the 4- vector part of the longitudinal polarization, we have an 0(e) 
suppression because the components of e° are of 0(m/E). The same suppression holds for 
the entire amplitude due to the selection rule (13.21) discussed above. In conclusion, these 
amplitudes are suppressed by a factor e. 

• Case (00 + ...+): This amplitude is even under (13.21) . and moreover the amplitude involving 
only positive polarizations and two scalar lines vanishes in the massless discussed 
in Section [2j We deduce that t is even and > 2 in (13. ip . To convince ourselves that it is 
actually t = 2, one can notice that the proof in the massless case relied on a specific choice of 
the reference vector, namely the momentum of a scalar line, and that there was no freedom 
to change it apart from choosing the momentum of the other scalar line. Analogously to the 
case ( — h ...+) discussed above, one can then see that these amplitudes must scale as e 2 . 
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e 3 , ...e n = +, ...+ e 3 , e 4 , ...e n = 0, +, ...+ 
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+ 
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+ 
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e 2 
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e 






e° 






e 







e 







e° 








e 2 








e 



Table 1: Scaling of the tree-level polarized amplitudes (ei,...e n ), in the axial gauge, with the pa- 
rameter e = m/E <C 1, see text. All the momenta are ingoing. The same suppression holds for the 
amplitudes related to the ones above by (+ O — ) in all the external states. 

• Case {ipip* + ...+): Going through the proof that these amplitudes vanish at tree level in the 
massless case, we see that it is valid also for a spontaneously broken gauge theory provided 
that the fermion is massless and we choose a light-like reference vector. If the fermion has 
a mass m^, we expect these amplitudes to be suppressed by only one power of m^/E. 

These findings are summarized in Table [TJ Notice that this choice of gauge is very convenient 
for tree level computations since the amplitudes (+...+), (0 + ...+) and {jpip* + ...+) all vanish 
(the last one only if = 0). Moreover each polarized amplitude is given by a sum of terms that 
are individually not parametrically larger than the total, i.e. there are never large cancellations 
among terms that individually grow faster than the total when the energy increases. 

4 Polarized amplitudes in covariant gauge 

In order to make contact with the more widely used covariant (Landau) gauge with d^A^ = 0, that 
implies e(p) ■ p = 0, we notice that the simplifications that occur in the axial gauge are crucially 
due to the structure (12. lft - (12. 31) of the polarization vectors. To understand the pattern of m/E 
suppressions in the polarized amplitudes in covariant gauge, it is then convenient to pass through 
a basis that resembles the axial one as much as possible. 
Recall that, for a generic on-shell 4- momentum: 

p M = (E p = \/m 2 + p 2 , p), (4.1) 

the polarization vectors corresponding to the three helicity eigenstates can be written as: 

«£(p) = -4(0, i P ±ij P ) , el(p) = -(\p\,E p p), (4.2) 

V Z 771 

where p = p/\p\ and i p , j p are two versors chosen in such a way that (i p , j p , p) form a right-handed 
ortho normal basis of the space-like 1R 3 . Of course, i p and j p are only defined up to a rotation along 
the p axis, and a different choice of the versors leads to the appearance of an overall multiplicative 
phase in e ± (p). 

As already said, in the context of massive vector bosons, helicity means "spin along the p 
direction" and it is well-defined when p ^ although it is not a Lorentz- invariant concept. It is 
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always implied that the helicities we consider are defined in a given reference frame, and we are 
particularly interested in the case of the CM frame. In fact the above definitions (14. 2 p are not 
covariant, and a generic Lorentz transformation mixes the different helicities. The helicity states 
are however not mixed by rotations, that is: 

R;e a (p) u <xe a (Rp)^ a = ±, (4.3) 

where R is a generic rotation, Rp is the rotated momentum, and we use the symbol oc to indicate 
the possible appearance of phase factors (see also Appendix [SJ. As for boosts, any boost along 
the p axis leaves e°(p) invariant, and it can either interchange or leave invariant the e ± (p) helicity 
states, depending on whether the sign of p is reversed or not. Finally, a boost not aligned with p 
mixes all the three helicity states. 

Now, let k be a generic 4-momentum such that p ■ k ^ and k ^ p; we want to construct two 
linear superpositions of the helicity eigenstates (14.21) that satisfy the conditions (12. 10 - 02. 3p . Since 
we are free to rotate the versors i p and j p: we can choose them such that j p ■ k — 0. Then it is 
easy to check that the 4- vectors that we look for can be written as: 

-+/ n 1 ( \/{k-p) 2 - {k\ + k 2 )m 2 \ , 1 k±m , , 

e±(p k) =- — y y k == — ± 1 ejAp) + _ ; e°(p) + 

2^ y/(k ■ p) 2 - k 2 m 2 J V2^/(k-p) 2 -k 2 m 2 

1 / ^J(k-p) 2 -(kl + k 2 )m 2 \ ^ (4 ' 4) 
2^ y/(k ■ p) 2 - k 2 m 2 T ) e ^ Ph 

where k± is the norm of the tranverse component of k with respect to p: 



k 



p-k) , fe = 4. (4.5) 



Notice that the second term in the right-hand side of eq. (14.41) becomes proportional to p M in the 
limit m — >■ 0, and that one of the two other terms vanishes. In this limit ( I4.4p simply expresses 
the fact that the polarization vectors ^(p, k) and e ± (p) represent the same physical state. The 
"axial-covariant" polarizations denoted by ^(p, k) can be thus defined in the standard covariant 
gauges as linearly- independent superpositions of the standard polarization vectors e ±,0 (p). As 
such, they can be completed to a full set of polarizations by defining a third vector orthogonal to 
the first two and to p, that we denote by e°(p, k). 

The crucial observation is now that a given polarized amplitude is gauge-invariant and it must 
depend only on the external states, that are specified by the definition of the polarization vectors. 
In particular, the relations (12. ip - (12. 31) are enough to fully specify which external states correspond 
to what we denote as "+" and "— " polarizations, and thus also the longitudinal one is fixed. As 
a consequence, a polarized amplitude computed in covariant gauge with the polarizations defined 
as in (14.41) must coincide with the corresponding amplitude computed in the axial gauge with the 
same reference vector. Notice that this equality is nontrivial since the Feynman rules are different 
and so is e°(p,k) ^ e°(p,k), as the first one it is orthogonal to p, the second one to k. Since 
e°(p,k) scales as E/m at high energy, we expect that large cancellations among subamplitudes 
are now present at least in the polarized amplitudes involving longitudinal polarizations, although 
the final result is the same due to the fact that gauge invariance is broken only spontaneously. 
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Table 2: Scaling of the tree-level polarized WW — > WW amplitudes with the parameter e — mjE, 
for £« 1, see text. The Higgs mass run is of order m. The missing combinations can be obtained by 
exploiting the C and P symmetry of the W Lagrangian. The polarizations used in this computation 
are the e(p,k) defined in (|4.4[) ; using instead the basis (|4.2[) one finds that the polarized amplitudes 
that here vanish become of 0(e 3 ) or 0(e 4 ), depending on whether a longitudinal polarization is present 
or not. See also |21) . 

To be specific, let us check these considerations in the phenomenologically interesting case of 
W — W scattering^]. We compute the scattering amplitudes using FeynArts/FormCalc [22] [23]. 
and we show in Table |5] the results for the various polarized amplitudes for the 2 — > 2 scattering 
using the polarization vectors ( 14. 4p . Taking into account that — > e T when an incoming external 
leg becomes outgoing^] and vice-versa, we see that the pattern of e l suppressions is exactly the 
same as the one that was deduced in the axial gauge in Table [TJ 

The only exception is (0, +; +, 0) ~ e 2 , while (0, +; 0, +) ~ e°. This is due to the conservation of 
the electric charge, that forbids some of the diagrams. Generally speaking, the suppression pattern 
of Tabled] is to be regarded as the "minimum amount of suppression", and it is always possible to 
have additional suppressions because of other reasons. In this case it is easy to understand why 
this amplitude is suppressed, using the methods of Sections [2] and |3] In fact in the configuration 
(0, +; +, 0) there is a change of electric charge in the scalar line, and then the only diagrams that 
contribute at tree level are those with two s — s — v vertices (where s means scalar and v means 
vector), and no v — v — v nor s — s — v — v ones. As a consequence this amplitude vanishes 
if we employ the momentum of one of the scalars as reference vector, as in the discussion of 
the (00* + ...+) case in Section |2] By using the arguments of Section [3] (the rule ( 13. 2 p and the 
continuity with the massless case), it is clear that the suppression becomes at least e 2 when we 
use a light-like reference vector. 

It is now interesting to ask which pattern we expect if we use the more conventional covariant 
polarizations (I4.2j) . This computation is performed in [21], and the result is similar to Table [2] 

9 We consider the SU(2)w interactions only, that is we set to zero the hypercharge and we do not include the 
QCD interactions 

10 In doing so one also changes the sign of the external momentum, so that an analytic continuation is understood 
in order to have positive energies. However one can also notice that the arguments of Sections [2] and [3] do not 
depend on the fact that the momenta are ingoing, and an e + can always be considered as an (e~)*. 
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with the only difference that the polarized amplitudes that here vanish become of 0(e 3 ) or 0(e 4 ), 
depending on whether a longitudinal polarization is present or not. This "double suppression" 
was noticed in [21], without explanation. It is clear that, on the contrary, passing through the 
axial gauge and our "axial-covariant" polarizations we have all the ingredients that are needed to 
explain this result. Consider in fact a polarized amplitude computed in covariant gauge using the 
polarization vectors e(p) defined in (14. 2p . The result has the general form: 

A {aM = A^ eL (4 P )* (4,oY , (4-6) 

where a,b,c,d = 0,±1. Take now this expression, that is linear in all the polarizations, and 
perform a change of basis by writing them in terms of the e(p, k). What we have to do is to invert 
the relations (I4.4p . but for the present purpose it is enough to notice that in the limit e C 1 
what happens is that the plus and minus polarizations mix at order e 2 , while the longitudinal 
polarization mixes with the transverse ones at order e. Notice that the high-energy behaviour of 
e° ,+ ' _ (p, k) is the same of that of e°' +, ~(p) respectively. As a consequence, the various polarized 
subamplitudes mix with each other, and more precisely A( a ,b;c,d) mixes with A( e ,f;g,h) at order: 

e |a-e|+|6-/|+|c-0|+|d-/i| j\ 

For example, A^^. equals zero, with the superscript indicating the basis of the e(p, k) with 

k 2 = 0, while A^^. + x is of 0(e 2 ). Since A^_ + . s in the basis of the e(p) takes contributions 

from *4.|^. H s at order e 2 , we conclude that it must be at least of 0(e 4 ) } as in fact it is since all 

the other contributions are of the same order or smaller. In the same way it is easy to see that 
all the "double suppressions" are understood, while the amplitudes that are "singly suppressed" 
(i.e. by e or e 2 ) remain of the same order in e. 



5 Example of phenomenological application: EWA 

The considerations of the previous Sections are of relatively general interest as they are valid 
in any theory with a gauge symmetry, unbroken (Section [2]) or spontaneously broken (Sections 
|3]and|4]). Let us now focus on the electroweak sector of the Standard Model, which is one of 
the most interesting subjects of study at present from the point of view of high-energy particle 
phenomenology. Indeed, uncovering the mechanism behind the spontaneous breaking of the elec- 
troweak gauge symmetry (EWSB) is arguably the main purpose of the CERN LHC. One of the 
most direct probes of the dynamics of EWSB is the high-energy scattering of the electroweak 
vector-bosons (W), especially the longitudinally-polarized ones. Since beams of W bosons do 
not exist, in order to study their interactions among themselves one has to rely on factorization, 
namely the fact that in a suitable kinematic regime the short- distance WW interactions can be 
separated from the well-known interaction of the W with the parton that emits it. As already said 
in the Introduction the resulting approximation, known as the effective W approximation (EWA) 
[20j[21j, is obsolete from the computational point of view since the experimental collaborations 
have enough computing power to simulate the exact process including radiative corrections. Nev- 
ertheless the EWA can still be a useful selection tool to understand which kinematic region is 
most sensitive to the dynamics of EWSB. 
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One can consider for instance the general process qX — > q'Y where q, q' are quarks and X, Y 
are respectively unspecified initial one-particle state and final state. Refering to [21] for details, 
we are interested in studying the interactions between q and X mediated by the W bosons and we 
speak directly of a quark since we have already factored-out the QCD processes^ that produce 
the quark from a hadron. Factorization in this case relies on the existence of a large separation 
of scales between the virtuality V 2 = m 2 — p\y of the collinear W emission and the hardness E of 
the relevant subprocess. In a practical situation, this amounts to require forward jets and large 
transverse momentum of the W bosons in the final state Y . Denoting the momenta of q and q' 
as: 



P q = (E, 0, 0, E) , P q , = (^(l-xyE'>+pi,p ± cos0,p ± sm ( p,E(l-x)) (5.1) 

the EWA amounts to say that in the limit m/i? C 1 and p± / E <C 1 the cross section for qX — > q'Y 
integrated over can be written as: 

where a denotes the polarization, Q stands for the charge, and Cq is a constant that depends 
on Q. The f a {x,p±) are computable splitting functions that describe the collinear emission of 
the equivalent W a whose momentum, that should be P q — P q ', is approximated to be on-shell in 
daiW^X^Y). 

Refering to [21] for discussions and generalizations, let us focus on the case in which X is an 
incoming W boson while Y consists on two W's. The proof of (I5.2p is based on a power expansion 
of the amplitude of the WW — > WW subprocess for small V 2 ~ m.&x{m 2 , p 2 ^} <C E 2 : 

A{V 2 )=A(0) + V 2 A'{0) + O(V 4 ) (5.3) 

and on the assumption: 

A'(0) ~ A(0)/E 2 . (5.4) 

This subamplitude represents the hard part of the process, the soft part being the collinear W 
emission from the quark that determines the splitting functions f a (%,P±) 111 (15. 2p . Concerning 
(15.41) . two important observations are in order. First of all one has to notice that, while .4.(0) 
is a gauge invariant quantity, the off-shell amplitude A(V 2 ) and its derivative -4/(0) are not. 
The reason is that, in order to obtain a gauge- invariant result, it is necessary to include the 
contributions from the "radiation diagrams" in which one of the final W's is emitted directly from 
the quark line. As a consequence, *4'(0) can be anything and actually in covariant gauges the 
condition (15 .4p is badly violated, apparently putting in doubt the validity of the EWA. In fact, 
although the diagrams that are neglected are not enhanced by any nearly-on-shell propagator, it 
can happen that in a "non-physical" gauge they give results that are much larger than the total, 
with large cancellations in the sum. The important point however is that, to prove the validity of 
the EWA, it is enough to find one gauge in which (15. 4p is valid, and it is shown in [2T] that the 
axial gauge does the job. The EWA is of course true in any gauge, but in some cases its proof can 
be more subtle. 



11 Meaning that the cross section that we write for an incoming quark has to be convoluted with the Parton 
Distribution Function of that quark in the colliding hadron. 
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Another important issue with (I5.3P is that, as discussed in the previous Sections, there are 
cases in which the on-shell amplitude -4(0) is suppressed by powers of e = m/E in some helicity 
configurations. Since V 2 is typically of order of the squared vector-boson mass m 2 , it may seem 
that the EWA (15.21) can fail in some cases for the polarized process. In fact if the on-shell amplitude 
-4(0) times splitting function were suppressed by more powers of e than the full amplitude, then 
one would have to conclude that the dominant contribution comes from the second term in (15. 3p . 
The EWA seems thus to be violated unless one is able to show that the full amplitude is always 
as suppressed as the on-shell WW amplitude times splitting function. 

Here is an example of a situation in which the machinery of Sections [2j|3] can be of use. In 
fact, besides a clear and simple understanding of the pattern of the e-suppression of the helicity 
amplitudes involving only vector bosons, we showed how one can easily understand also the 
suppression of the amplitude involving two fermions and gauge bosons with same helicity. To 
prove the validity of the EWA in a channel in which there is some amount of e-suppression for all 
the polarizations of the equivalent W, such as q(W)~ —> q' (W) + (W) + , it is then sufficient to show 
that the full amplitude, including the femions, is suppressed by the same power of e — m/E as 
the dominant WW — >■ WW subamplitude. Since this suppression is independent of V 2 , it affects 
the entire A(V 2 ) in 15.31 and thus also *4'(0). 

Consider the above example of q{W)~ — > q'(W) + (W) + with massless quarks and the corre- 
sponding effective-iy process (W e ff) a {W)~ —> (W) + (W) + . From the previous Sections we know 
that, in the axial gauge with light-like reference vector, the amplitude for the full process van- 
ishes. On the other hand, looking at the subamplitude times splitting amplitude, when a = — 
the subamplitude vanishes while when a = 0, + it is the splitting amplitude that can be made to 
vanish by suitably choosing the reference vector as the 4-momentum of one of the two fermions. 
Notice in fact that, considering for example the splitting amplitude with a left-handed fermion 
t/>_(Pi) /. a (k,p)ip_(p2), one can make the case a = — vanish by choosing p — pi, or alternatively 
the case a = + vanish by choosing p = p 2 (see (I2.15P ). not both simultaneously. The amplitude 
with a = instead vanishes in both cases. 

The situation becomes more involved if one employs the usual covariant-gauge definition of 
the polarization vectors, denoted by e(p) in Section HI In this case we know that the scaling of 
the amplitude A a of the process (W eff ) a (W)- ->• (W) + (W) + is: 

A+ ~ e 2 , A ~ e 3 , A~ ~ e 4 (5.5) 

and in the splitting amplitudes there is an additional e suppression in the case a = 0. The validity 
of (15. 2p is not obvious now, unless one is able to show that the full tree-level amplitude including 
the fermions is suppressed at high energy at least by a factor e 2 . Without such a proof, in [2T] it 
was checked numerically that in this and similar situations the EWA remeins nevertheless valid, as 
always up to corrections of order max{m 2 , p 2 ± } / E 2 . Using our results it is immediate to understand 
why the EWA holds: the full amplitude for q(W)~ — > q'(W) + (W) + is zero in axial gauge, while 
it is unsuppressed if we flip the helicity of one of the W's. On the other hand we know by (14.41) 
that the polarizations e ± {p) in covariant gauge have a component of e T (p,k) of order e 2 . As a 
consequence, the full amplitude is suppressed by e 2 for any value of the virtuality, and it matches 
the behaviour of the EWA result. To be precise, one may question that in covariant gauge eq. 
(15.41) does not hold in general, so that the second term in (15.31) could be of the same order as the 
first one. However it can be shown that, by suitably modifying the gauge, one can make (15 .4p valid 
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also in a certain class of covariant gaugeo. Similar situations can be understood analogously, 
sometimes making use of considerations like those made in in Section H] about the case in which 
some diagrams are not present because of the conservation of the electric charge. 

6 Conclusions 

In this paper we studied some properties of the helicity amplitudes for spontaneously broken gauge 
theories at high energy The concept of helicity becomes frame-dependent in the case of massive 
gauge bosons of mass m, and in general it can be "violated" by a larger amount with respect to the 
massless case. Our main result is the generalization to the massive case of the vanishing theorems 
that are valid in the case of unbroken gauge theories. In particular we provided a set of "selection 
rules" that specify by which power of m/E a given helicity amplitude is suppressed at energies of 
order E 3> m. Moreover we did it using only elementary means, and we also presented a simple 
discussion of the known results that are valid in the unbroken case. To show the usefulness of 
our approach, we discussed how one can employ our methods to understand some aspects of the 
effective W approximation in the polarized case. 
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A Axial gauge 

Refering to [20] and references therein for more details, the axial gauge is defined by the gauge- 
fixing condition: 



that is enforced by a delta function in the functional integration. An important feature of this 
choice is that, unlike in the covariant gauges, the mixings among Goldstones and gauge fields 
do not cancel out. As a consequence the propagator is nondiagonal when regarded as a matrix in 
the 5-dimensional space (A^, s) spanned by the four components of plus the Goldstone boson 

12 This alternative proof of the EWA will be presented in (24) . 
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s. Its explicit form is: 



Vu{q,k) = — — , (A.2) 

\f ( n b) - -a + q ^ kv + ^ - g a 



Af ss (q,k) 



{qkf 
m 2 k 2 



(qkf 



By construction, the projector matrix Mu annihilates the vector (fc M , 0) and has thus rank 4 for 
general q. At the pole q 2 = m 2 , the rank is further reduced down to 3, and the projector can be 
rewritten as: 

Afjj(q 2 = m 2 ,k)= J2 ^(q,k)e x j(q,kr (A.3) 

A=+,-,0 

where the polarization vectors ef(p, k) have vanishing fifth component and satisfy fl2.1|) - fl2.3p while, 
for A = 0, the 4- vector €j(p, k) describes the longitudinal polarization: 

el(p,k) = (A.4) 



(pfc) 

e° s (p,k) = - 



m?k 2 



(pk) 2 

Let us consider the transformation properties of these polarization vectors under rotations. In 
order to be more concrete, we adopt a "constructive" approach. First of all, from the definition 
(IA.4j) it is clear that the longitudinal one transforms as: 

el(Rp,Rk) = Rtl{p,k) (A.5) 

where R is any rotation matrix. For the transverse states we are particularly interested in rotations 
along the axes specified by p when p 2 = (for any k 2 ), and rotations along k with k 2 = (for any 
p 2 ). In general however eq. (14. 3 p holds. 

Let us start with a rotation R p e of an angle 9 along the particle momentum p. Without loss 
of generality, we can take p to be along the third axis and moreover, since R p e commutes with a 
boost along the same axis, we can go to the frame in which k has no third component. The two 
4-vectors transform then as: 

R% = Rl(p o ,0,0,\p\)=P, (A.6) 
R? e k = Rg(k , ki, k 2 , 0) = (k , ki cos8 + k 2 sin 9, k 2 cos 9 — ki sin^, 0) . (A. 7) 

To construct e ± {p 1 k) we need to choose two 4-vectors a and b with a 2 = b 2 = —1/2 such that 
{a, b, k} form an orthogonal basis of the 3-dimensional subspace orthogonal to p, and such that 
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(a x b) ■ p > to ensure right-handedness. In our case we can define I (for p = 0, i.e. po = \p\): 

■f^^ 1 '^' • m =T2 ( l'°' lk i ] - (A - 8) 

Up to an overall phase factor we can define the transverse polarization vectors satisfying fl2.ip - fl2.3p 
as e ± (p, k) = a(k) ± ib(k). It is then immediate to verify that: 

e ± (Rlp, R$k) = a(Rlk) ± ib{R p e k) = e ±ie R p d (a(k) ± ib(k)) = e^R^fa k) [p 2 = 0] . (A.9) 

Consider now the case of a rotation along the direction of the reference vector k with k 2 = 0. 
Everything proceeds as above, with the roles of p and k interchanged, except for the choice 
of a and b. In fact, with the definition (lA.8p . one has (a x b) = (— ki/ko, — fe/^o, 1) so that 



(a(k) x b(k)) ■ p > but when the construction is done with k p one has to choose the opposite 
ordering for a{p) and b{p). As a consequence: 

e ± (Rjp, Rjk) = b(Rjp) ± ia{R % e p) = e Tie Rj(b(p) ± ia(p)) = e^R^ip, k) [k 2 = 0] . (A.10) 

To consider the more general cases in which we rotate for example along p with p 2 ^ 0, one 
can notice that the expression for a(k) and b(k) in (IA.8P can be generalized to: 



ru\ 1 f kl fi _l x\ i _l x k * x klk2 hP°fi i ^\ 



(A.ll) 



where: 



^w)"'- 1 - (A ' 12) 

One can then see that the same transformation rules are still valid. 

To summarize, the relevant transformation properties of the transverse polarization vectors 
can be written as: 

Jjf^Cp,*) = e^e ± (Rtp = p,R{k), (A.13) 
R^fak) = e ±ie e ± {Rjp,Rjk = k) . (A.14) 

Notice that the fact that the rotation along the k axis brings a phase which is opposite to that 
from a rotation along the p axis is what one expects from (12. 4p and (12. 5p . 
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17 



[2] 

[3] 
[4] 
[5] 

[6] 
[7] 

[8] 
[9] 
[10] 

[11] 
[12] 

[13] 
[14] 
[15] 
[16] 
[17] 
[18] 
[19] 
[20] 
[21] 
[22] 
[23] 
[24] 



Z. Kunszt, Nucl. Phys. B 271 (1986) 333. 

M. T. Grisaru, H. N. Pendleton and P. van Nieuwenhuizen, Phys. Rev. D 15 (1977) 996. 
M. T. Grisaru and H. N. Pendleton, Nucl. Phys. B 124 (1977) 81. 
S. J. Parke and T. R. Taylor, Phys. Rev. Lett. 56 (1986) 2459. 
F. A. Berends and W. T. Giele, Nucl. Phys. B 306 (1988) 759. 

Z. Bern, L. J. Dixon, D. C. Dunbar and D. A. Kosower, Nucl. Phys. B 425 (1994) 217 

| arXiv:hep-ph/94 03226] 



E. Witten, Commun. Math. Phys. 252 (2004) 189 arXiv:hep-th/03 '12171 



R. Britto, F. Cachazo and B. Feng, Nucl. Phys. B 715 (2005) 499 arXiv:0412308 



R. Britto, F. Cachazo, B. Feng and E. Witten, Phys. Rev. Lett. 94 (2005) 181602 \airXw :0501052. 
N. Arkani-Hamed, F. Cachazo, C. Cheung and J. Kaplan, JHEP 1003 (2010) 020 \arXiv:0907.5418\ 



N. Arkani-Hamed, J. L. Bourjaily, F. Cachazo, S. Caron-Huot and J. Trnka, JHEP 1101 (2011) 
041 \arXiv:1008.2958\ 



M. L. Mangano and S. J. Parke, Phys. Rept. 200 (1991) 301 arXiv:hep-th/0509223 



L. J. Dixon, In *Boulder 1995, QCD and beyond* 539-582 arXiv:hep-ph/9601359 



L. J. Dixon, J. M. Henn, J. Plefka and T. Schuster, JHEP 1101 (2011) 035 \arXiv: 1010. 39911 



M. E Peskin, arXiv.1101.2414 



R. H. Boels and C. Schwinn, Phys. Rev. D 84 (2011) 065006 arXiv: 1104.2280 



S. Dittmaier, Phys. Rev. D 59 (1998) 016007 arXiv:hep-ph/ '9805 445 



V. V. Andreev, Phys. Rev. D 62 (2000) 014029 arXiv:hep-ph/0101140 



Z. Kunszt and D. E. Soper, Nucl. Phys. B 296 (1988) 253. 



P. Borel, R. Franceschini, R. Rattazzi and A. Wulzer, JHEP 1206 (2012) 122 arXiv: 1202. 1904 



T. Hahn, Comput. Phys. Commun. 140, 418 (2001) arXiv:hep-ph/0012260 



T. Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118 (1999) 153 \arXiv:hep-ph/ 9807565 
A. Wulzer et al, unpublished, in preparation. 



18 



